gravity data inversion using l1-norm stabilizer

Authors

کیلان راست بین

دانشگاه آزاد همدان سعید وطن خواه

دانشگاه تهران وحید ابراهیم زاده اردستانی

استاد گروه فیزیک زمین

abstract

in this paper the inversion of gravity data using l1–norm stabilizer is considered. the inversion is an important step in the interpretation of data. in gravity data inversion, the goal is to estimate density and geometry of the unknown subsurface model from a set of known observation measured on the surface. commonly, rectangular prisms are used to model the subsurface under the survey area. the unknown density contrasts within each prism are the parameters which should be estimated. the inversion of gravity data is an example of underdetermined and ill-posed problem, i.e. the solution can be non-unique and unstable. thus, in order to find an acceptable solution regularization should be imposed. solution is usually obtained by minimizing a global objective function consisting of two terms, data misfit and the regularization term. data misfit measures how well an obtained model can reproduce the observed data. usually, it is assumed noise in gravity data is gaussian, therefore a l2–norm measure of the error between observed and predicted data is well suited for data misfit. there are several choices for a stabilizer, depends on type of features one wants to see from inverted model. a typical choice is a l2 –norm of a low-order differential operator applied to the model, which also a priori information and depth weighting can be incorporated (li and oldenburg, 1996). in this case the objective function is quadratic, then minimization of the function results a linear system to be solved. however, the models recovered in this way are characterized by smooth feature which are not always consistent with the real geological structures. there are situations in which the sources are localized and separated by sharp, distinct interfaces. to deal with this problem, during last decades, researchers have proposed a few types of stabilizer. last and kubik (1983) presented a compactness criterion for gravity inversion that seeks to minimize the area (or volume in 3d) of the causative body. portniaguine and zhdanov (1999) based on this stabilizer, who named the minimum support (ms), developed the minimum gradient support (mgs) stabilizer. for both constraint, the regularization term can be written as the weighted l2–type norm of the model. therefore, the problem of the minimization of the objective function can be treated same as conventional tikhonov functional. the only difference is that a priori variable weighting matrix for model parameters incorporated in the regularization term. thus the iteratively reweighted least square (irls) algorithm is required to solve the problem. other possibility for stabilizer is the minimization of the l1-norm of model or gradient of model, the latter indicates total variation regularization. the l1–norm stabilizer allows occurrence of large elements in the inverted model among mostly small values. therefore, it can be used to obtain sharp boundaries and blocky features. although the l1–norm stabilizer has favorable properties, in reconstruction of sparse models, its numerical implementation in a minimization problem can be difficult because its derivatives with respect to an element is not defined at zero. to overcome this difficulty, in this paper, the l1–norm stabilizer is approximated by a reweighted l2 –norm term. the algorithm is extended to gravity inverse problem, which needs depth weighting and other priori information to be included in the objective function. for estimating the regularization parameter, which balances between two terms of objective function, the unbiased predictive risk estimator (upre) method is used. the solution of the resulting objective functional is found using generalized singular value decomposition (gsvd), also provides for efficient determination of the regularization parameter at each iteration. simulation using synthetic data of a dipping dike demonstrates that the method is capable to reconstruct focused image, boundaries and slop of the reconstructed model are close to those of the original model. the method is applied on gravity data acquired over the gotvand dam site, in the south-west of iran. the results show rather good agreement with those obtained from the boreholes.

Upgrade to premium to download articles

Sign up to access the full text

Already have an account?login

similar resources

Tomographic inversion using l1-norm regularization of wavelet coefficients

We propose the use of l1 regularization in a wavelet basis for the solution of linearized seismic tomography problems Am = d, allowing for the possibility of sharp discontinuities superimposed on a smoothly varying background. An iterative method is used to find a sparse solution m that contains no more fine-scale structure than is necessary to fit the data d to within its assigned errors. keyw...

full text

An alternative transformation in ranking using l1-norm in data envelopment analysis

Jahanshahloo et al. (Appl Math Comput 153:215–224, 2004) propose a method for ranking extremely efficient decision making units (DMUs) in data envelopment analysis (DEA) using super-efficiency technique and l1-norm and they show that the presented method is able to eliminate the existing difficulties in some methods. This paper suggests an alternative transformation to convert the nonlinear mod...

full text

A method for 2-dimensional inversion of gravity data

Applying 2D algorithms for inverting the potential field data is more useful and efficient than their 3D counterparts, whenever the geologic situation permits. This is because the computation time is less and modeling the subsurface is easier. In this paper we present a 2D inversion algorithm for interpreting gravity data by employing a set of constraints including minimum distance, smoothness,...

full text

3D gravity data-space inversion with sparseness and bound constraints

One of the most remarkable basis of the gravity data inversion is the recognition of sharp boundaries between an ore body and its host rocks during the interpretation step. Therefore, in this work, it is attempted to develop an inversion approach to determine a 3D density distribution that produces a given gravity anomaly. The subsurface model consists of a 3D rectangular prisms of known sizes ...

full text

Non-linear stochastic inversion of 2D gravity data using evolution strategy (ES)

In the current work, a 2D non-linear inverse problem of gravity data is solved using the evolution strategies (ES) to find the thickness of a sedimentary layer in a deep-water situation where a thick sedimentary layer usually exists. Such problems are widely encountered in the early stages of petroleum explorations where potential field data are used to find an initial estimate of the basin geo...

full text

My Resources

Save resource for easier access later


Journal title:
فیزیک زمین و فضا

جلد ۴۲، شماره ۲، صفحات ۳۳۷-۳۴۸

Keywords
in this paper the inversion of gravity data using l1–norm stabilizer is considered. the inversion is an important step in the interpretation of data. in gravity data inversion the goal is to estimate density and geometry of the unknown subsurface model from a set of known observation measured on the surface. commonly rectangular prisms are used to model the subsurface under the survey area. the unknown density contrasts within each prism are the parameters which should be estimated. the inversion of gravity data is an example of underdetermined and ill posed problem i.e. the solution can be non unique and unstable. thus in order to find an acceptable solution regularization should be imposed. solution is usually obtained by minimizing a global objective function consisting of two terms data misfit and the regularization term. data misfit measures how well an obtained model can reproduce the observed data. usually it is assumed noise in gravity data is gaussian therefore a l2–norm measure of the error between observed and predicted data is well suited for data misfit. there are several choices for a stabilizer depends on type of features one wants to see from inverted model. a typical choice is a l2 –norm of a low order differential operator applied to the model which also a priori information and depth weighting can be incorporated (li and oldenburg 1996). in this case the objective function is quadratic then minimization of the function results a linear system to be solved. however the models recovered in this way are characterized by smooth feature which are not always consistent with the real geological structures. there are situations in which the sources are localized and separated by sharp distinct interfaces. to deal with this problem during last decades researchers have proposed a few types of stabilizer. last and kubik (1983) presented a compactness criterion for gravity inversion that seeks to minimize the area (or volume in 3d) of the causative body. portniaguine and zhdanov (1999) based on this stabilizer who named the minimum support (ms) developed the minimum gradient support (mgs) stabilizer. for both constraint the regularization term can be written as the weighted l2–type norm of the model. therefore the problem of the minimization of the objective function can be treated same as conventional tikhonov functional. the only difference is that a priori variable weighting matrix for model parameters incorporated in the regularization term. thus the iteratively reweighted least square (irls) algorithm is required to solve the problem. other possibility for stabilizer is the minimization of the l1 norm of model or gradient of model the latter indicates total variation regularization. the l1–norm stabilizer allows occurrence of large elements in the inverted model among mostly small values. therefore it can be used to obtain sharp boundaries and blocky features. although the l1–norm stabilizer has favorable properties in reconstruction of sparse models its numerical implementation in a minimization problem can be difficult because its derivatives with respect to an element is not defined at zero. to overcome this difficulty in this paper the l1–norm stabilizer is approximated by a reweighted l2 –norm term. the algorithm is extended to gravity inverse problem which needs depth weighting and other priori information to be included in the objective function. for estimating the regularization parameter which balances between two terms of objective function the unbiased predictive risk estimator (upre) method is used. the solution of the resulting objective functional is found using generalized singular value decomposition (gsvd) also provides for efficient determination of the regularization parameter at each iteration. simulation using synthetic data of a dipping dike demonstrates that the method is capable to reconstruct focused image boundaries and slop of the reconstructed model are close to those of the original model. the method is applied on gravity data acquired over the gotvand dam site in the south west of iran. the results show rather good agreement with those obtained from the boreholes.

Hosted on Doprax cloud platform doprax.com

copyright © 2015-2023